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Abstract 

The paper proposes an algorithm for regularization of the self-energy expressions for a 
Dirac particle that meets the relativistic and gauge invariance requirements. 

Within the second order formulations of the "old" perturbation theory for free motion, 
the expression for the upper integration limit q is a slowly enough varying function of 
particle impulse. For a particle at rest, q = m; for p/m in the range from to 1, 




4 6 \ 

.00864-^-r - 0.00217-^r ; 

to 4 77i° / 



for the ultra-relativistic case |p| S> m : q « 1.5m. 

For 4D perturbation theory, on introduction of the limiting 4-impulse, L 2 = Lg — L 2 , it 
is shown that with a large time component, Lq/tu 3> 1, the spatial values of L» are limited 
and are the same as the components of the introduced limits of integration qc L 2 = q 2 . 

Within the proposed algorithm, in the second-order of the perturbation theory, the 
renormalized Dirac particle mass is 

m (2) = m + Ato (2) = to ^1 + 1.115— \ , 

where too is the bare mass of the particle. 

The electromagnetic radius of particle at rest determined by equating the self- energy 
and the electrostatic energy e 2 /r em is about half the Compton wavelength, r em w l/2m. 



1 




It is commonly known that the calculation of some physical effects results in undefined 
diverging expressions, when the perturbation theory formalism of existing formulations of field 
and particle quantum theories is used. If the quantum-field theories are renormalizable, the 
undefined expressions are removed in all orders of the perturbation theory by renormalizing 
mass and charge of relevant particles. In particular, quantum electrodynamics is such a theory. 

In the lifetime of the relativistic quantum-field theories, many researchers tried to solve 
the regularization problem of infinite expressions appearing in the computations from various 
viewpoints. The problem, however, remains unsolved. 

This paper suggests, as its authors believe, a natural way for solving the problem of finite 
renormalization of particle mass by the example of quantum elecrodynamics. 

First, use the formulations from the "old" perturbation theory given by Dirac for quantum 
electrodynamics in p]. 

The self-energy operator in the second order of the perturbation theory is: 

n An J ^\E\ + Ikl - vE* Ik| 

In (pj) and hereinafter fi = c = 1; 
a\ (3 are Dirac matrices; p = p + k; 



E = an + (3m + eA ; n l = p l - eA l ; n % = p % + k l - eA l \ 

E = an + (3m + eA - p = (m 2 + p 2 ) 1/2 ; p = (m 2 + p 2 + 2pk + k 2 ) V2 ; 

A , A 1 are scalar and vector potentials of the external electromagnetic field. 
Expression (JTJ implies the averaging over v : (1/2) (Y(v = 1) + Y(y = —1)). 
In action on state \i) with energy Ei, 

E \i) = Ei \i) , Y \i) = Y \i) (2) 

Dirac considered three modes for the regularization of expression (JTJ in [T]. 
Mode a . Modulus of k is limited, |k| < q, where q is some number. 

Mode b . Here the sum of energy moduli of three particles involved in the interaction event 
is limited, \E\ + |k| + \E^\ < 2q. 

Mode c . Here the sum of kinetic energies of three particles is limited, 

\E - eA Q \ + |k| + \Ei - eA Q \ < 2q. 

For large q, the correspondence is established: 



Y a = Y b + \an + \eA ; 
o 2 



Y a = Y c + -an 



(3) 
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As a result, Dirac prefers mode c as a technique allowing us to eventually obtain relativistically- 
and gauge-invariant expressions for the self-energy and vacuum polarization. As this is always 
the case in the relativistic quantum field theories, the numerical value of q is not determined. 

This paper uses mode a to demonstrate the possibility to obtain the relativistically- and 
gauge-invariant expression for the self-energy with simultaneous estimation of the upper inte- 
gration limit q. 

First, consider the case where there are no external electromagnetic fields (A = 0, A 1 = 0). 
Expression (|T} becomes 



Yn 



Air* 



a, 



up + ap + (3m dk 



M Pobo + |k| - vpo) |iv. 

Consider a non-relativistic charged particle having impulse |p| Cm. 
On summation over /x and integration over azimuthal angle (p, we obtain 



(4) 



<i i 



Yn 



e 

7i 



-1 



-upo + a(p + k) + 2(5m 
Po(Po + |k| - vpo) 



kd(cos 9)dk. 



(5) 



Expanding the integrand in © over the particle impulse |p| up to quadratic expressions 
p 2 /m 2 and integrating give: 



Yn 



2e 2 



7T 



m 



ln(rr) 



+ 



m 



ki(x) 



+ (3m 



In (a;) + 



3(x 2 

x 2 



1 m 2 p / 

x 2 2m 2 \3 x 2 + 1 
1 m 2 \ 



m 



4 m 2 



+ 



Am 2 



ln(x) 



1 m T 
24^2" 



m 



2m 2 \ 3x 2 + l (x 2 + I) 2 3 (x 2 + 1) 



x = k + Vm 2 + fc 2 , A ■ m = q + Jm 2 + q 2 . 



A-m 



Use equalities 



— m — — — h ap + (3m ] \i) 
2m I 



P 2 P 2 

(i\ ap \i) = (i\ — \i) ~ (i\ — \i) 
p m 



to obtain 



(i\Yo\i) 



e 

7T 



(3m 



2 ln(A) + 4^"4 

e! I 5 4 

+ m I 3{A 2 + 1) + (A 2 + l) 2 



3(A 2 + 1) 3 + e)}'^ 



(6) 
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Reasoning from the relativistic invariance condition for the self-energy operator Y , require 
that the coefficient of p 2 /m in @ vanish. The natural condition for the upper limit of inte- 
gration q in (JHJ follows herefrom: 

5 4 8 1 n 

+ 7TS 7T77 - 77^7 + 77 = Or 



3(A 2 + 1) (A 2 + l) 2 3(A 2 + 1) 3 6 
9 ( Q 2 1 



(g 2 + m 2 )V2 V q 2 + m 2 2 



0. (7) 



The solution to equation (j2J) is g = ±m, in addition to the trivial one, q = 0. 

The value g ~ m means that the upper integration limit imposes the restriction on the 
distances. They should be longer than Compton wavelength of Dirac charged particle: |x| > 
1/m, which is quite reasonable from the viewpoint of quantum mechanics. 

The value q ~ m for the upper integration limit is also reasonable in terms of "Zitterbewe- 
gung" effect for Dirac particle. In fact, if Dirac Hamiltonian is written in the form 

H D = an + [3m + eA ; n = p — eA 

in the presence of external fields, then the operator of coordinate x % can be represented as 
follows 0: 



x i = xi + 7T*— — — - — -t - \aie~ 2i{HD ~ eAo)t - 1 



H D -eA 4 u (H D -eA c 



where a l is the matrix integration constant. 

The third term of x l operator highly oscillates with frequency ~ 2m. This term can be 
represented as 

Ax' = - ( a 1 - 7T € — 1 , ) TT 1 , (8) 
2 V H D -eAj H D ~eA v ' 

It is seen from © that Ax* ~ 1/m, and from this standpoint it is quite natural to restrict 
the consideration to distances |x| > 1/m. 

Note that when previously estimating the self-energy of the non-relativistic charged particle 
in the self-field electrodynamics in Foldy-Wouthuysen representation, one of the authors also 
came to the necessity of the finite upper integration limit with A « m 

Clear, the value q ~ m obtained from the expansion of © to quadratic values p 2 /m 2 can 
vary when the following expansion terms are included. The value of q can also change with 
external electromagnetic fields A^(x) introduced. 

Consider the function 

e 2 /3 1 1\ 

(i\ F(p, q, m) \i) = (i\ Y \i) - — (i\p m [- ln(A ) + - - ] \i) , (9) 
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where A = A\ q=m = 1 + ^/2. 

The second addend in © is relativistically invariant. In its form, the addend corresponds 
to the infinite renormalization term to be added to the particle mass in standard quantum 
electrodynamics, if the integration limit q — ► oo in A(q). 
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Table 1: The upper integration limit q(p) 



p/m 





0.01 


0.02 


0.04 


0.05 


0.06 


q/m 


1 


1.000000000088 


1.0000000014 


1.000000023 


1.000000055 


1.00000011 


p/m 


0.07 


0.08 


0.09 


0.1 


0.2 


0.5 


q/m 


1.00000021 


1.00000036 


1.00000058 


1.00000088 


1.000014 


1.000512 


p/m 


0.8 


1 


10 


100 


500 


1000 


q/m 


1.00297 


1.00647 


1.26716 


1.448 


1.48117 


1.494 



The function F(p, m, q) is the relativistically non-invariant expression that (with our ap- 
proach) should vanish with an appropriate choice of the upper integration limit q. 
Table [T] gives q/m versus p/m, for which F{p,m,q) = 0. 
The following equalities were used in the numerical computations: 

<*l fa 10 = <«1 ( (m2 + p2) i/ 2 ) 10 > <*i «p io = <*i ( K+ P p2) i/ 2 ) 10 • 

The computed data shows that the upper integration limit q is a slowly enough varying 
function of p/m. Within the range of p/m values from to 1, the function is amenable to 
approximation by the polynomial expression 

q » m (l + 0.00864^- - 0.00217^-r ) . (10) 
V m 4 m 6 J J 

For the ultra-relativistic case of |p| ^> m, the function approaches q ~ 1.5m asymptote. 
Thus, when choosing the upper integration limit q(p) according to Table [2 for free motion, 
the average value of the self-energy operator (J3]) is a relativistic invariant, 

(i\ Y |i) = - (i\ (3m (\ ln(Ao) + - -\ \i) = 1.115- (i\ (3m \i) . (11) 

We now turn to the modern 4D relativistically invariant formulation of the perturbation 
theory in the quantum field theory. 

The self-energy operator (or mass operator) in the second-order of the perturbation theory 

is 

^Tci 2 f 2m — p + k 4 



M(p) = -- — -e 2 / — — -- — ■ — —r—d k. (12) 

yFJ (2tt) 4 J [{p - k) 2 - m 2 }k 2 K ' 

For the external electron line of Feynman diagram, 

p 2 = pi — p 2 = m 2 , p\i) = m \i) . 

Introduce, according to the notion of the limiting 4-impulse, L 2 = L\ — L 2 , and perform 
the integration in (|T2T) over the variable k° using the Feynman rule of pole bypass. For L /m ^> 
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1 the integration result is 



M(p) 



An 2 



P / dk 





k 


— ak + f3m 




k|(p + |k 


-Po)(po ~ 


k 


+ Po) 





k + ak — (3m 


' |k| fa- 
te 


k 


-Po)(Po + l k l + Po) 
p ) — «k + /5m 


Po(Po + |k 


-Po)(Po - 


k 


-Po) 



(p + Po) - ak + (3m 
Po(Po + |k| +Po)(Po - |k| +p ) 



It can be shown that on the algebraic transformations the average of the expression in state \i) 
equals to that of expression (pjj for Y : 

(i\M(p)\i) = (i\Y \i). 

Hence, the spatial components of the introduced limiting impulse L are equal, Lj = qi. Since 
for free motion the only relativistic invariant is particle mass m, then L 2 = L\ — q 2 = C ■ m 2 , 
where C is a numerical factor. As Lo/m 3> 1, the factor C should be also much higher than 
one. With q varying as a function of the particle impulse, the time component L should also 
vary accordingly, with the invariant L 2 remaining invariable. 

Note that if the 4D integration in (|T2jl is performed in the conventional manner, then 



expression (|12j) can be written as 

/ x e 2 
M(p) = — m 

7T 

At the limit L' 2 > m 2 




m 1 
arctan 



L' 2 /m 2 



U 81 + L' 2 /m 2 



(13) 



M(p) 



e 

f« — m 

7T 



/2 1 
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In 



+ 



Ln 2 +L' 2 ) because 



which coincides with a similar expression obtained in [4]. 

The components of impulse V in (fT3|) form the Euclidean 4D sphere [L' 2 — ^ 
of the rotation of the path of integration over variable k° in the complex plane. Besides, the 
components of impulse L' differ from those of the previously introduced impulse L because of 
the change in the domain of integration from the replacement of variables, k — p — > k 0]. 

To determine L', once again use equality 



Then, 



/2 s 



In 1 + 



L 



(i\M(p)\i) = (i\Y (q)\i). 

1 L' 2 /m 2 



L' m 

arctan 

2m U 



81 + L' 2 /m 2 
4/3, L' 2 w 1. 



1,115. 



(14) 



The solution to equation (fT4]) is L'/m = 1.335 ~ 4/a, u " ^ i.sm' 
With consideration of the above-described regularization of the expression for the self-energy 

of free Dirac particle, the particle mass change due to the renormalization in the second order 

of the perturbation theory is 



Am' 2 ' = — m 

7T 



^ln(A )+ ' 



AA 2 



1.115 — m. 

7T 



(15) 
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The re-normalized particle mass is 



m 



m + Am (2) = m 1 + 1.115 




(16) 



where m is the bare mass of the particle. 

By setting the self energy of Dirac particle at rest, (i|Yo(|p| = 0) |z), equal to electro- 
static energy, e 2 /r em , we obtain the electromagnetic radius estimate for non-relativistic particle: 
r em « l/2m. 

Given external fields A M (x), the integration limit q and time component L will depend on 
their magnitude with retained invariant L 2 . 

For weak external fields v4^(x), where in self-energy expression (JTJ we can restrict ourselves 
to inclusion only of terms not higher than quadratic in the generalized particle impulse, the 
results of the calculation for the anomalous magnetic moment of the particle and Lamb shift 
of energy levels in the second order of the perturbation theory within the approach developed 
in this paper will be presented in the following publication. 

In the presence of high external electromagnetic fields, electromagnetic radius r em will reduce 
owing to changes of the integration limit q, thus ensuring quantum electrodynamics applica- 
bility up to ultra-short distances. 

Of course, it would be desirable that ^(^(x)) be obtained from first principles in the future 

Summarize the results of the paper. 

1. An algorithm for regularization of self-energy expressions for Dirac particle satisfying the 
requirements of the relativistic and gauge invariance is proposed. 

2. The expression for the upper integration limit q is a relatively slowly varying function 
of particle impulse within the second-order formulas of the "old" perturbation theory for free 
motion. 

For a particle at rest, q = m; for p/m in the range from to 1, 



for the ultra-relativistic case |p| ^> m, q ta 1.5m. 

It is shown for 4D perturbation theory upon introduction of limiting 4-impulse L 2 = Lq — L 2 
that with a large magnitude of the time component, L /m ^> 1, the spatial values of Lj are 
limited and are the same as the components of the introduced limits of integration qi. L 2 = q 2 . 

For the conventional integration techniques for Feynman integrals in the second order of the 
perturbation theory the limiting integration impulse is evaluated: L' 2 = L' 2 + L' 2 ps 16m 2 /9. 

3. In the second order of the perturbation theory, the re-normalized mass of Dirac particle 



qmm[l + 0.00864^— - 0.00217 



m 




is 



m 



m + Am (2) = m 1 + 1.115 




where m is the bare mass of the particle. 
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4. The electromagnetic radius of particle at rest that can be found by equating the self- 
energy and the electrostatic energy, e 2 /r em , is approximately half the Compton wavelength, 
r em m I /2m. 
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